QCD traveling waves beyond leading logarithms 
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We derive the asymptotic traveling-wave solutions of the nonlinear 1-dimensional Balitsky- 
Kovchegov QCD equation for rapidity evolution in momentum-space, with 1-loop running coupling 
constant and equipped with the Balitsky-Kovchegov-Kuraev-Lipatov kernel at next-to-leading loga- 
rithmic accuracy, conveniently regularized by different resummation schemes. Traveling waves allow 
to define "universality classes" of asymptotic solutions, i.e. independent of initial conditions and of 
the nonlinear damping. A dependence on the resummation scheme remains, which is analyzed in 
terms of geometric scaling properties. 



I. INTRODUCTION 

o 

o . 

In the large- AT C and "mean-field" approximation of high energy (high density) QCD, the density of gluons with 
transverse momenta k in a target evolves with rapidity Y according to the nonlinear Balitsky-Kovchegov (BK) equa- 
, tion. This equation is supposed to capture essential features of saturation effects. In the 1-dimensional approximation, 
^"*^ ■ and within leading logarithmic (LL) accuracy, it reads P, 0, 0] 

dyJV(L, Y) =a X (~d L ) Af(L, Y) - aAf 2 (L, Y), (1) 

(N . 

where L — log(fc 2 /fcg), with being an arbitrary constant. In the LL approximation the characteristic function of 
• the kernel has the standard Balitsky-Kovchegov-Kuraev-Lipatov (BFKL) form [HID, Q, namely 

Jo ' X(7) = 2^(1) - #y) - iKl - 7), (2) 

m ■ 

f^*) ' and the coupling a = a s N c /ir is kept fixed. 

The goal of the present paper is to extend the known traveling-wave method @, H, H[ for obtaining asymptotic 
solutions of the nonlinear equation ((T|) to the case where one considers the equation at the next-to-leading logarithmic 
| (NLL) accuracy. In that case, the QCD coupling constant is running, i.e. 

^ m = ~, b= ^-2N f 

' . V ' bL 1 \2N C ' V ' 

Ph. 

^ ' and the 1-dimensional BK equation which we consider reads 



, — i 

X 



bLd Y Af(L,Y) = x(-d L ,d Y ,a)Af(L,Y) - Af 2 (L,Y) . (4) 

As we shall recall later on, x(~ 9l, dy, a) is a "renormalization-group improved " NLL kernel. It follows from specific 
resummation scheme implying higher order contributions while keeping the known expression [13, HH of the NLL 
term in the kernel. Note that there could be also NLL corrections to the nonlinear term, but, as we shall see from 
the traveling-wave properties, they do not change the asymptotic solutions. 

To summarize the present situation of the QCD traveling waves for our purpose, it has been shown @, H, that the 
BK equation with fixed coupling a and in the LL approximation belongs to the same universality class as the Fisher 
and Kolmogorov-Petrovsky-Piscounov (F-KPP) equation [T3, EH- This means in particular that the BK equation 
admits solutions in the form of traveling waves M{L — v g aY). L has the interpretation of a space variable while 
t — aY , interpreted as time, is an increasing function of rapidity Y. v g is the critical velocity of the wave, defined 
in this case as the minimum of the phase velocity. From the point of view of QCD the traveling-wave solution for 
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the quantity Af translates into the property of geometric scaling [3] which yields N (k 2 / 'Q 2 S (Y)) , where the function 
Qs(Y) cx exp (vgaY) is called the saturation scale. Hence the asymptotic traveling-wave solutions of the BK equation 
satisfy geometric scaling (up to sub-dominant scaling violations which we will not analyze in the present work). 

The existence of this solution depends on the general form of the initial condition Af(L,Yo) for large values of L. 
It is however universal, i.e. independent on the details of the kernel or the specific form of the nonlinear term (e.g. 
independent of NLL nonlinear terms) provided that it plays the role of a damping saturation term. This unique 
condition is that N(L,Yq) decreases at large L at least as rapidly as an exponential exp(— 70K) with 70 > j c , where 
7 C is the critical anomalous dimension, solution of the equation 

X(7c) = JcX'ilc) ■ (5) 

The corresponding value for -f c obtained with the LL kernel is j c = 0.6275 while the value for the critical velocity 
v g = x'ilc) = 4.883. The above condition is fulfilled for high energy QCD due to the "color transparency" property 
of the gluon density i.e. Af(L, Y — Y ) ~ e~ L , valid at large L. Hence, 70 = 1 > "f c ■ 

Within the framework of the traveling- wave approach the results for the gluon density Af(L, Y) and the saturation 
scale Q 2 (Y) were obtained in the LL approximation Moreover, in ||, the method was also extended to the 

case of the equation with the LL kernel and 1-loop QCD running coupling constant ([3|). In fact one considered the 
equation similar to (U|), but keeping the LL kernel ©. In this case, the traveling- wave method leads to a different 
universality class of solutions than the F-KPP one, with, in particular, geometric scaling property in y/Y instead of Y. 

One would like to extend the above procedure the for BFKL kernel considered at the NLL accuracy and find the 
analytic solution in the framework of the traveling- wave approach for the corresponding BK equation This has 
been done for the case of the fixed coupling a [151 ] . The case of the running coupling, which has not been considered 
so far, is the subject of this paper. 

As we know now, although the NLL corrections were calculated for the BFKL kernel [l(| [H| they turned out to 
be negative and so large that they are of no use for a s if it is not extremely small. However, several equivalent 
ways to cure this pathological behavior were proposed. They are based on the observation that the problem of the 
NLL corrections comes from the existence of spurious collinear singularities. This singularities may be canceled by 
resummation of the collinear terms at all orders satisfying, at the same time, the renormalization group constraints. 

In practice, various resummation schemes were developed from which we will use the S3 and S4 schemes [16j and 
CCS scheme [TtJ • It is important for further use (and possible generalization to other schemes), to distinguish two 
types of resummation among the kernels. In the CCS "implicit" scheme [l7|, the higher order resummation appears 
through the dependence of the kernel on two variables only x(l = ~ u> = dy), where uj — 0{q) drives the higher-order 
corrections. In the S3 and S4 schemes [l6[, an "explicit" dependence on the coupling constant appears, leading to a 
triple- variable dependence xil — ^ — dy , a) . As we shall see later on, this introduces significant analytical (and 
eventually phenomenological) differences in the traveling wave solutions. 

The aim of this paper is to apply the traveling-waves method with running coupling at the NLL level. As an 
outcome we obtain the result for the gluon density and the saturation scale valid at large Y which is universal and 
can be used with any resummed NLL kernel. The specific cases of the S3, S4 and CCS kernels are studied in more 
detail. We also compare the traveling-wave solutions with that obtained with a method of the linear BFKL evolution 
in the presence of absorbing boundary conditions [l8|], which has been applied using the CCS scheme (l9j |. 

The plan of the paper is the following. In Sec. HI] we present in detail the calculation which leads to the solution of 
the BK equation with the NLL BFKL kernel in the case of running coupling. We arrive at the analytic asymptotic 
expressions for the saturation scale and the gluon density. This allows us to define in which universality class the 
traveling-wave solutions lie, depending on the type of the resummation scheme. In Sec. IIIII we specifically analyze 
the dependence of the logarithmic derivatives of the saturation scale on the resummation scheme used. We check the 
consistency and generality of our approach by comparison with the previously known results for the CCS scheme. 
Finally, the conclusions and outlook are given in Sec. IIVI 

II. BK EQUATION WITH NLL BFKL KERNEL AND RUNNING COUPLING 

Following the general method [8j we first write the solution to the linearized version of the BK equation at NLL ([¥]). 
It has the form of double Mellin transform [ItJ 

m > Y) = J £ / £ M,(7 ' u) exp (~ 7L +ujY+ ~L x{i > u) ) ' (6) 



with 



X( 7 ,w)= [dr/xW,u) , (7) 

and 7 being an unspecified constant. Indeed, using the saddle-point method for the integration over 7 at large enough 
L, one gets the saddle-point equation 

- L +T~ X(7.") = 0, (8) 

or equivalently in operator form 

bLd Y N{L,Y) = X (-d L ,d Y )M(L ) Y) , (9) 

which is nothing else than the restriction of the BK equation Q to its linear part. 

Here again \ as we U as X may depend on a. In order to solve Eqs. ([6]) and (O, we shall make the approximation of 
fixing a in the kernel at some, phenomenologically motivated, value, knowing that it applies only to the higher-order 
corrections. We shall come back to this approximation later on. In order to simplify the notation we do not write the 
dependence on a explicitly. Let us remind that the CCS kernel, as well as all other "implicit" schemes do not include 
this dependence on the coupling constant. 

As a next step we perform the saddle point integration over lu in the limit of large Y and obtain 

M(L,Y) = J |j-A/- (7) ex P (- 7 £ + F(u; s )r), (10) 

where F(lu s ) — Y b u, (2^(7,^) — u> s X(^f,u> s )^ and the condition for the saddle point uj s is given by the implicit 
equation 

Ybuj 2 s ~ X{ 1 ,uj s )+uj s X{ 1i oj s ) = . (11) 

We introduce the notation in which the "dot" means the derivative with respect to to whereas the "prime" means the 
derivative with respect to 7. Let us now expand the integral of the kernel ([7]) near lu = 

x (l ,„) = ££!!^. (12) 

p=0 p ' 

Using a similar expansion for its derivative X with respect to to and substituting both quantities into Eq. (|lip we 
have 



F6+ix( 7 ,0) 



^ = X (7,0) -|E^h)T X(P) (7.0)^| . (13) 



by collecting terms by powers of co s . It is clear from Eq. (fl3"|) that for asymptotic Y, tu s ~ y -1 / 2 while in the 
convergence domain of the series the remaining terms between braces are at most of order Y~^l 2 . As we shall check 
later on, they may contribute only to non-universal sub-asymptotic terms. We also checked this by looking for an 
iterative solution where the kernel x(7j w ) is expanded around loq and truncated at the order P. The universal terms 
are shown not to depend on either loq or the truncation. 

A. Traveling-wave critical parameters 

Considering then (|13|) up to the order two, we obtain 



*<*°> (14) 



r6+iX( 7 ,0) 



The saddle point value behaves like ui s ~ Y 2 . With this form of u> s the gluon density is given by 

Af(L, Y)= j ^M(l) exp(- 7 L + n( 7 )t), (15) 
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FIG. 1: Graphical determination of the critical exponent y c for the three resummed NLL kernels. The curve corresponding to 
the CCS scheme (at u = 0) coincides with the LL curve. 



where time is interpreted as 



t = VY + Yo, (16) 



with Yq = X( r y 1 0)/2b. Interestingly, this non-universal term in formula (|14p absorbs the arbitrary constant 7 in 
X(7, 0). The dispersion relation has the form 



O( 7 ) = ^X( 7 ,0). (17) 

In analogy to the LL case, following 8], a critical group velocity (defined as the minimum of the phase velocity in the 
wave language) is obtained as 

v g = n{ lc )/ lc = n'( 7c ) . (is) 

However, 7c determined in such a way still depends on the arbitrary constant 7. Thus, requiring v g to be independent 
on the choice of 7 means imposing du g ( 7 )/c? 7 = = du ff ( 7c )/d 7c . This is because the dependence of the velocity on 
7 comes through 7c only. Applying this condition to Eq. I|17p one gets 



dv g (j)/d^ = d(n( lc )/ 7c )/d Jc =0 v g = y ^ 7 7 c e ' 0) , (is) 

eliminating all dependence on the arbitrariness in the definition of X in Eq. 0. This is to be contrasted with the 
situation for the linear problem with NLL kernels. 

The value of j c at the NLL level in general depends on the resummation scheme. This is clearly shown in graphical 
form in Fig. [TJ Geometrically, the value of "f c is given by the tangent to the characteristic function of the kernel, 
different for each NLL scheme. Note also that the curve corresponding to the CCS scheme at u> = is nothing else 
than the LL curve given by Eq. ([2]), and thus the critical parameters are the same in this case. This holds for any 
other "implicit" scheme which recovers the LL kernel at ui = like e.g. those proposed in [2(| HH ■ 



B. Asymptotic solution of the BK equation 



The linearized version of the BK equation @ with the kernel expanded around loq up to the second order and the 
rapidity variable Y changed to time variable t = \JY + Y is given by 

"LL" 



h J- dtM = {-^s i + i x "(a£ + 2 7c a i + 7c 2 ) 



kw-kxM-y ^ + w^-\ 9i) k (2o) 



"NLL" 

We singled out twoparts in the above equation. The part denoted as "LL" has already been present in the LL case 
(cf. Eq. (33) from Q). The remaining part called "NLL" contains new terms which originate from the dependence of 
the NLL BFKL kernel on uj. In accordance with the approach developed in Q we take the ansatz for the solution in 
the form [HJ 

Af(z > t)=t a G(z)exp(- lc zt a ), z= L-v g t + c(t) ^ (21) 

where we assume c(t) — /3t k ^ 1 and the constants a, (3 and k are to be determined. In order to recover the LL equation 
in the limit x,x\x ~> 0, we set a — | and k = ^, as in [8]. In fact this ansatz remains valid beyond the leading 
order. It is easy to check the consistency of this choice by looking at the time dependence of the new terms generated 
at higher order in Eq. ([20)1 . When all derivatives which appear in "NLL" part are calculated it turns out that their 
leading contributions are proportional to £3. This means however that the terms in "NLL" part contribute only at 
the order t~$ since each derivative is multiplied by at least the factor t~ l . Thus, the BK NLL linearized equation in 
this approach has exactly the same form as in the LL case and reduces to the Airy equation 

G"(*) = ^^(*-4/3)G(*). (22) 

The condition G(z) ~zasz->0 allows to fix the constant j3 to 



„ Sl , (23) 



4 V lcV g b 

where £i = —2.338 is the zero of the Airy function. Finally the result for the gluon density is given by 



and the saturation scale up to a multiplicative constant has the form 



Here x{lc, 0) is the NLL BFKL kernel resummed in a given scheme (and possibly taken at some value of a, see Sec.H]). 
Hence, the solution of the BK equation with the resummed NLL kernel and running coupling has the same functional 
form as the solution for the LL kernel given in [8j- In particular, in the saturation scale the leading exponential term 
proportional to the time variable t is supplemented by a universal term in t*, sub- leading by order t~s . Note that the 
higher order contributions to the saddle-point Eq. (|12p. being sub- leading by the order Y~? ~ t _1 , do not interfere 
with this analysis and are expected to be non-universal. 

Despite this formal similarity with the LL case with running coupling, the critical exponents are not the same for 
LL and NLL case if x(7,w = 0, a) ^ x(l)i given in ||5J). This is due to the fact that the NLL result depends on the 
resummation scheme used for the kernel xili w = 0, a). However, the traveling- wave solutions are not sensitive to the 
values of the kernels for u> ^ 0, contrary to the applications of the NLL BFKL linear equation without saturation 
(see, e.g. 0,11). 




III. SATURATION SCALES BEYOND LEADING ORDER 



As we have seen, the time variable (JT^J) scales like t ~ Y 1 / 2 . For a quantitative analysis of the solution, one defines 
the logarithmic time-derivative of the saturation scale (or intercept) as 

M ( Y) - dios(q - ( ;; )/A ^ ) . (26 , 

In order to compare with the usual definition of geometric scaling implying a linear rapidity dependence of the 
saturation scale, one may also consider an "effective" intercept 

dlog(Q2(y)M2 ) 

X ff(Y) = BWs \ " , (27) 



dY 



so the relation between them is A'(Y~) = 2t \ e Jf(Y). 
From Eq. (f2T)|) we obtain 



A . (F) X'U.0) &H| (28) 



where the time £ is defined by Eq. p^|) . In Fig. [2] we show A* from Eq. (|2"5)) as a function of F 1 / 2 for the three different 
resummation schemes S3, S4 and CCS and the value of the coupling a = 0.15. The result depends on the scheme 
used. For Y 1 / 2 — > oo the logarithmic derivative A* goes to its asymptotic value equal to the group velocity v g , which 
is also scheme-dependent. 



A. CCS scheme 



In fact, there appears a difference between the scheme CCS and the schemes S3 and S4. In the first case, since 
the higher order effects in the kernel are all given as a function of u>, one has identically x(7,o> = 0) = x(w), i.e. the 
LL kernel of In mathematical terms, this means that the CCS scheme falls into the same universality class of 
solutions as the equation with a LL kernel and running coupling constant. 

In order to check this result, we compared the traveling- waves solution with that obtained with the different method 
applied to the CCS scheme [ijj]. In this method, one considers only the linear evolution term supplemented by 
the absorbing boundary conditions [l8| . It has been confirmed in @, H, [§| that the direct traveling- waves approach 
applied with the LL kernel leads to the same asymptotics (the traveling- wave method allows to add the third universal 
term Q for the fixed coupling case while an eventual third term for the running case remains an open problem). 
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FIG. 3: Comparison of the logarithmic derivative Xl obtained in the traveling-waves approach with the result of [IS 
sake of compatibility we use CCS scheme and Yo = 0. The corresponding formulas are given in Eqs. ([28} and (|29[) . 



For the 



Extracting the asymptotic analytic form of the solution found in 19] (cf. Eqs. (58) and (59) therein) , and after 
changing it to our notations, one obtains 



'2x(7c,0) 



b^fc 




X"(7c,0) 
v /2 7c 6 X ( 7c ,0) / 



6 t~ 



2b 



7cX(7c,0) 



(29) 



We see that the two results (|28|l and ([2"9"| are consistent up to the corrections of the order t~a . However, the corrections 
at this order are not expected to be universal. The detailed comparison of the two results for A* given by Eqs. 
and (|29|) is shown in Fig. (3[ We observe that they converge for the asymptotic values of Y. 
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FIG. 4: Comparison of the logarithmic derivative A s obtained in the traveling-waves approach in the CCS scheme and various 
Yo values with the results of [13]. The notations BFKL "L in uj" , "NL in w" correspond to the running coupling with the LL 
kernel and two different treatments of the CCS scheme, respectively, see [l!| for details. 



For a quantitative comparison at lower rapidities, we can also confront our result for \ e /f with that given in p^ ]. 
This is presented in Fig. H) We see that \ e J* based on the definition of Eq. (f2"T|) with appropriate Yg can successfully 
mimic the results from 19]. In our analysis, varying Y$ plays the role of parameterizing typical non-universal terms, 
i.e. terms which depend on the initial conditions, details of the kernel, or of the method used for extracting the 



asymptotic behavior. This is a hint for the compatibility of the results of [191 ] with the universality class defined by 
the LL kernel with running coupling constant. The analysis of non-universal terms is beyond the scope of the present 
paper and deserves further study, since these terms may be phenomenologically important. 



B. S3 and S4 schemes 
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FIG. 5: Dependence of the logarithmic derivative obtained in the traveling-waves approach in the S3 and S4 schemes as 
a function of a. Top: S3 scheme; Bottom: S4 scheme. The curves go up when the value of a decreases. The fixed curve 
corresponding to the CCS scheme, i.e. the same as for LL kernel, is also indicated in both cases. 



The schemes S3 and S4, by contrast, are not expected to lie in the same universality class as the previous one. 
Indeed, the resummed kernels depend explicitly on the the value of the coupling constant. This is depicted in Fig. [SJ 
where one can see how the time-derivative of the saturation scale varies with a. In order to characterize the universality 
class of these schemes one would have to go one step beyond the approximation of a fixed a in the kernel that we 
used in the beginning. In fact one would have to first solve exactly the linear operator equation 

bL d Y N{L, Y) = X (-d L , dy,a(L) = 1/bL) Af(L, Y) , (30) 

and then apply the general traveling-waves method to find the critical velocity and asymptotic solution. We postpone 
this detailed study for future work. What, however, is indicated by our study is that the "universality class" of the 
solution may be different from the CCS case. It remains to be found which kind of geometric scaling property will be 
satisfied. 

Hence we obtain the new result from the QCD traveling-waves approach that the specific asymptotic solutions of 
the BK equation at NLL accuracy depend on the resummation schemes. In particular they vary between the schemes 



which include an explicit or implicit dependence on the running coupling constant in the definition of the NLL kernel. 
This introduces a theoretical distinction between NLL effects. It would be an interesting issue to know whether this 
distinction may appear at physically reachable rapidities. This deserves certainly more study in the future. 

IV. CONCLUSIONS 

We considered the 1-dimensional Balitsky-Kovchegov (BK) equation in momentum space with running coupling 
constant. The Balitsky-Kovchegov-Kuraev-Lipatov (BFKL) kernel was taken at the next-to-leading logarithmic (NLL) 
accuracy with higher orders following known resummation schemes and satisfying renormalization-group (RG) con- 
straints. 

Using mathematical properties of nonlinear equations, we derived the traveling-wave solutions, which are valid in 
the limit of large rapidities Y and obey universality properties, i.e. are independent on the specific form of the initial 
conditions, the detailed form of the kernel and the nonlinearities. 

We found that the results for the gluon density and the saturation scale acquire dependence on the resummation 
scheme. For a scheme with an implicit higher order dependence, i.e. where the NLL kernel is a function of 7 = — 8l 
and uj = 8y only, one is expected to fall into the same universality class as in the LL, running coupling case. This is 
explicitly derived for the CCS scheme [l~7l |. 

For schemes with an explicit higher order dependence, such as S3 and S4 of [la ] . where in addition the NLL kernel is 
an explicit function of a, the asymptotic solution does change as a function of the value of a. The precise determination 
of the universality class when a is also considered as running in the NLL kernel is an interesting challenge for the 
future work. 

It is interesting to address further questions about the properties of QCD saturation beyond leading logarithms. For 
instance, it will be useful to ask which is the universality class of the BK equation (or whether they are modified) when 
different ways of implemen ting the running coupling constant are used, since there remain an ambiguity concerning 
this problem (see e.g. [Ill Hyj). On a more phenomenological ground, one would like to say something about the 
ph ysic ally reachable rapidity region where non-universal terms may be important. For instance, a way was proposed 
in [27l . |28| to take into account the actual form of nonlinearities beyond the asymptotics where their specific form does 
not play a big role (except in the selection of the critical velocity by setting the unitarity limit) . Finally, there remains 
the problem of a more complete QCD solution with NLL accuracy going beyond the mean-field approximation. We 
postpone these studies for further work. 
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